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ABSTRACT 

The comiection between renormalons and power corrections is investigated for 
the typical infrared renormalon integral assuming the effective coupling constant 
has an infrared fixed point of an entirely perturbative origin. It is shown the full 
answer differs from the Borel sum by a power correction. 
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It is generally believed that the large order behavior associated to infrared (IR) 
renormalons, which makes QCD perturbation theory "non Borel summable", reflects 
an inconsistency related to the Landau ghost. In this note, I examine how the renor- 
malon problem is resolved when the coupling constant is free of Landau singularity 
and approaches a non trivial IR fixed point at small momenta. I assume the fixed 
point arises entirely within a perturbative framework, through higher order pertur- 
bative corrections. I will show that in this case, IR renormalons are also present, but 
the exact result differs from the Borel sum by a (complex) power correction, which 
removes all inconsistencies. 

Let us first review the standard argument [1,2] for renormalons. Consider the 
typical IR renormalon integral : 




aef}{k/Q,a) (1) 

where a is the coupling at scale Q in some arbitrary renormalization scheme, and 
aeff{k) a renormalization group (RG) invariant effective coupling (I assume n > 0, 
so that the integral in eq.(l) is IR convergent order by order in perturbation theory). 
Let us compute the perturbative Borel transform R{z). If R{a) has the formal power 
series expansion : 



p=0 



R{z) is defined by : 



Ri-) = E ^ (2) 

The series eq.(2) are believed to have a finite convergence radius (at the difference of 
those for R[a)), and therefore allow for an all order definition of R{z). One can then 
define an all order "perturbative" resummed Rpt{oc) by the Borel representation : 



RpT{a) = dz exp (^-^^ R{z) (3) 
It is convenient to express R{z) in term of the Borel transform of cte// : 



aeff{k/Q,a) = dz exp (^-^^ aeff{k/Q, z) (4) 

Inserting eq.(4) into cq.(l), and interchanging the order of integrations, one indeed 
recovers eq.(3) with R{a) = Rpj ioi) and : 



R{z)=l n—i^-j a,ffik/Q,z) (5) 

IR renormalons arise as an IR divergence [1] of the integral in eq.(5), resulting from 
the IR behavior of aeff{k/Q, z) (which follows solely from the RG invariance of CKg//) : 



/9: 
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aeff{k/Q,z) ~ aeff{z)exp -z y3o In 

e < Q2 (g) 
where (3^ and f3i are (minus) the one and two loop beta function coefficients. Using 



eq.(6) into eq.(5) one indeed finds for z 



n/(3o : 



R{z) ~ aeffiz ~ Zn f^^^^^ ^1 _ \)i+s (7) 

with S — j^Zni i.e. R{z) displays a cut singularity [2] at the renormalon position z — 
Zn, which generates according to eq.(3), since Zn > 0, an 0{exp{—Zn/(x)) imaginary 
part in RpT{a). In the standard case where cce// has a Landau ghost at some scale 

< Q^, this fact causes no surprise, since the defining integral eq.(l) itself involves 
an ambiguous integration over the Landau singularity. But a paradox arises if Oe// 
has an IR fixed point : then R{a) should be perfectly well defined, with no imaginary 
part according to eq.(l), whereas eq.(3), together with eq.(7), still yield an ambiguous, 
imaginary amplitude for RpT{a) ! But eq.(7) is correct in both cases : as mentionned 
above, it follows solely from RG invariance and the representation eq.(5) (which is 
presumably always valid since it is correct order by order in perturbation expansion in 
z, and the corresponding series are convergent). Thus IR renormalons are also present 
in the IR fixed point case. The only possible conclusion is that eq. (3) is not in fact a 
valid representation of eq.(l), i.e. that R{a) in eq.(l) does not coincide with Rpri^Oi) 
in eq.(3) in the latter case. This is despite the fact I assumed eq.(4) does instead 
correctly represents (at least for large enough k) aeff{k) (which is also assumed to 
have no renormalons), which can thus be determined in principle from "all order" 
perturbation theory, and qualifies the present framework as being "perturbative" . 

Indeed, there is a loophole in the previous derivation that R{a) = Rpt{q:), which 
is actually already present when there is a Landau singularity. Assume e.g. Oe// is 
the one loop coupling : 



aeff{k/Q,a) 
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Then [1] : aeff(k/Q, z) — exp[—z Po In (^)] 6q.(5) yields 



R{z) = / 
Jo 



n- 
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The problem with the previous argument is that eq.(4) itself is a valid representation 
of aeff{k) only if k is not too small. Indeed in the present example the integral in 
eq.(4) converges at 2; = oo only if ^ + /3o In (^-^^ > 0, i.e. for k^ > A^. For /c^ < A^, 
aeff{k) has to be represented instead by a Borel integral over the negative z axis : 



aeff{k/Q,a) ^ - [ dz exp (--\ a^ffik/Q, z) (8) 

Similar remarks apply in the general case, where eq.(6) implies, provided aeff{z) 
decreases no faster then exp{—cz) at large z, aeff{k) will have a Landau singularity 
at some scale k'^ — A^ (see however the fixed point case), below which the Borel 
representation eq.(4) will break down, and eq.(8) should be used instead. 

These observations suggest that the correct procedure in the Landau ghost case is 
to first take a < (and slightly complex if l3i 7^ 0), which implies < A^, so that in 
the whole integration range in eq.(l) one has k^ < < A^, and the representation 
eq.(8) is valid throughout the range (this means choosing a in the domain of attraction 
of the trivial IR fixed point). Manipulations analoguous to those performed above 
are now justified, and yield : 



R{a) = -j^ dz exp (^-^ R{z) (9) 

Finally, R{a) for a > (where > A^) is defined as the analytic continuation of 
eq.(9), which yields the standard Borel representation eq.(3) with R = Rpt- 

In the case where aeff{k) has an IR fixed point, a similar problem arises if it 
happens again that the representation eq.(4) is not valid below some scale k = k„iin, 
even if kmin does not correspond to a Landau singularity this time. This is possible 
if, as in the Landau ghost case, aeff{z) does not decrease too fast at large z (the 
latter assumption seems necessary, because a too fast decrease of aeff{z) (e.g. [1] if 
aeff{z) = exp{—cz^)), although insuring the validity of eq.(4) for all k's, usually yields 
an aeff{k) which blows up too fast as k 0, making the original integral eq.(l) IR 
divergent). Furthermore, and this is the crucial difference with the Landau ghost case, 

for k < krain ■> 

aeff{k), which remains positive, wiU not admit a Borel representation 
over the 2; < axis either. (Note that for /c — > 0, aeff{k) is always in the strong 
coupling region, since it approaches a non-trivial fixed point, at the difference of the 
Landau ghost case, where for k ^ one enters an (IR) asymptotically free region 
below the Landau singularity, since no other fixed point is available). The above 



derivation that R — Rpt can thus not be extended to the IR fixed point 
expected. Let us now give a general argument that these two functions actually differ 
by a power correction. Splitting the integration range in eq.(l) at /c = kmim one has : 



_^ , , / rmn 



(k \ fQ' dkUk^Y (k ^ 



= R_ + R+ 

Eq.(4) cannot be used inside the low momentum integral R-, which can however be 
parametrized as a power correction since : 

\ J k'^ \kminj \kmin J 

( kL„ \ " 



where cimm = o:{Q = kmin), and RG invariance has been used in the first step. On 
the other hand, using eq.(4) in R^, one gets : 



r r f z\ fQ' de fey 

R^ = J dzexp^-- J n — ^- 



To go further, it is necessary to know the k-dependence of aeff{k/Q, z). This can be 
done easily in the special case where f3i = 0, if one chooses a{Q) to be the one loop 
couphng. Then one can show[l] that : aeff{k/Q, z) = aeff{z)exp[—z(3o ln{k'^/Q'^)], 
and one gets : 



|2 ( k'^y^ 



k^ \Qy w -"1 
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Hence : 



z 



.2 \ " 



= Rpria) - J RpT^amin) 

One therefore ends up with the result : 



R{a) = RpT{a) + [-^\ [i?(a^.n) - i?pr(am^n)] (10) 



where the coefficient of the power correction is given by the discrepancy between the 
exact ampUtude and its Borel representation. Eq.(lO) is equivalent to the statement 
that R{a) = RpT{<y) + const / {Q'^)"' , and is also correct when f3i ^ (see below). 
Note it is not possible the power correction vanishes (as it does in the Landau ghost 
case) since R{amin) is real, whereas RpT{cimin) is complex due to the effect of the 
renormalon (the power correction must be complex to cancell the imaginary part of 
RpT{a)). 

I now illustrate the previous discussion with the example of the 2 loop coupling : 



dink?' 



(3,{a,ffY - (3,{a,^fY (11) 



I shall consider both the standard case i3i/(3q > where there is a Landau singularity, 
and the case (3i/f3o < where aeff{k) has an IR fixed point at am = —(3o/f3i (which 
actually occurs in QCD for a large enough number of flavors). Remarkably, R{z) 
can be computed exactly with a straightforward change of variable, adapted from a 
similar one suggested in ref.[3]. Defining the Borel variable by : 



1 



Zn l + |a 
with a = C(eff{k = Q), and using the solution of eq.(ll) : 



(12) 




^T- ^ h— + - + ^ (13) 




(together with the similar relation at k — Q) one obtains : 

dk'^ f k^Y exp{-^ - 

n-rir I I aeff{k) — —dz- " 



^2 yQ2j '-e/M'V (l-^)'""' 

which suggests the looked for Borel transform is : 



CXp I —ttZ. 



To complete the proof, it remains to determine the new integration bounds 



1) Assume first /3i//?o > : then a^ff has a Landau singularity at the scale 
A2 = K^exp (-^In (1^) + J)- I assume > A^, so that a > 0. For k'^ = Q^, 



cte// = a, and z — ; ai — , a^ff — oo, and z — zl = — < Zn ; decreasing 

k'^ below A^, z becomes complex ; finally, for A;^ 0, Oe// 0~, and z — > +oo. R{a) 
in eq.(l) then takes the form of a Borel integral along a path in the complex plane, 
which divides into two complex conjuguate branches at z = zl, and can be deformed 
to above or below the positive real axis to yield the standard Borel representation 
eq.(3) with R{z) as in eq.(14). 



2) On the other hand, if Pi/Po < 0, the situation is actually simpler : as k'^ 
decreases from to 0, Oe// increases from a to the IR fixed point am, and z 
increases from to Zn through real values so that eq.(l) becomes : 



(the renormalon singularity is integrablc at z = Zn, since S < now). We therefore 
check that R{a) is not given by the Borel sum RpT^a) in the IR fixed point case. 
Rather, it is given by Rpt^o) minus an exponentially small 0{exp{—Zn/cx j) correction : 



Z 
1+5 



\ Z„J \ ZnJ 

= RpT + Rnp (16) 

The latter is just a (complex) power correction, in accordance with the general 
expectation. Indeed, eq.(15)-(16) simplify by performing the change of coupling : 
1/a— 1/a + Pi/ Pq, and one gets in particular : 



C exp {—Zn/a){—l/aY 

c{-iY (AVg')" (17) 

with C = ^r(l — 5){znY , and cq.(13) at k = Q was used in the last step (in this ex- 
ample, the subtracted term —Rnp corresponds exactly to the "minimal" prescription 
of ref.[4] to "regularize" IR renormalons) . The same method can deal with the case 
a < 0, where one is in the domain of attraction of the trivial IR fixed point (provided 
the condition 1 + > is also satisfied if (3i/(3o < 0). Then a^ffik) monotonously 
increases from a to 0~ as k'^ decreases from to 0, hence z < and decreases from 
to -oo. Thus : 



R 



NP 



— [ dz exp ( — \ 

Jz„ ^ ^ (\ — 

V Zn 
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R{a) = f dz exp (--) '^^^^ (18) 



oo \ a 



i.e. a Borel integral over the negative z axis, in agreement with eq.(9) (note that 
eq.(15) is not the analytic continuation of eq.(18)). 

For completness I give also the result if aeff{k) satisfies the RG equation : 

daeff _ -/?0 "e// 



dlnk^ 1 - ^a^ff 

where the inverse beta function has only two terms. The appropriate change of 
variable in this case turns out to be precisely the one suggested in ref.[3] : 

z _^ a 



(with a = aeff{k = Q)), which yields for /3i//3o > (using also the convolution 
theorem) : 

R[a) — a H / dz exp I 



whereas for (3i/(3q < 0, where ccg// has an infinite IR fixed point, one gets 

5 1 
R{a) — a H / dz exp ( — - 



6'l + 6Jo'^\aJ {l-^y^' 



and shows that in this case too (with a different C) : R{a) — RpT{o()—C{—iy{A^/Q^)"'. 

For an arbitrary /3e// function with an IR fixed point one expects however the 
answer to be of the more general form : 



R(a) = RpT(a) + (AyQ^n-C(-lY + C) (19) 

where C and C are real, and independent of Q, and C is proportionnal to the renor- 
malon residue. It may be worth mentionning that taking the derivative of both 
sides of eq.(19), and going to Borel space, one easily derives a relation between R{z) 
and the Borel transform of (3eff{a) — da/ din (with a = CKeff{k — Q)) : 



R{z) - ^z R{z) - - r dy b{z - y)y R{y) 
n n Jo 



(20) 



where b{z) is the Borel transform of b{a) = _ Eq.(20) can be used to 

rederive the previous two results, and allows to deal with more complicated examples 
as well. Note also that C and C, being independent of Q, have dropped from eq.(20), 
and in fact any Q-dependence in these coefficients would be inconsistent with the 
assumption of the perturbative nature of the coupling (equivalent here to assuming 
that Peff^oi) coincides with its Borel sum). 
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